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Abstract. - We consider counterions in the presence of a single planar surface with a spatially 
inhomogeneous charge distribution using Monte-Carlo simulations and strong-coupling theory. 
For high surface charges, multivalent counterions, or pronounced substrate charge modulation 
the counterions are laterally correlated with the surface charges and their density profile deviates 
strongly from the limit of a smeared-out substrate charge distribution, in particular exhibiting 
a much increased laterally averaged density at the surface. 



The recently revived interest in charged soft-condensed matter systems reflects that there 
are still many open questions, despite the enormous amount of work available in this field. One 
example of high experimental relevance is the discrete nature of charged surface groups, or, 
more generally, the inhomogeneity of substrate charge distributions, and how it affects various 
thermodynamic properties in an aqueous environment such as forces between charged particles 
or the counterion distribution. The importance of the discreteness of charged surface groups 
has been experimentally established in colloidal fiocculation[l] and deposition studies [2] and 
was recently reviewed [3]. On a much larger length scale, chemically micropatterned substrates 
with charged and neutral patches can be used for controlled colloidal deposition [4] and DNA 
immobilization^]. Theoretically, charge- modulated surfaces have been studied by various 
mean-field approximations [6, 7, 8, 9, 10], liquid state theory[ll, 12], as well as computer 
simulations[13, 14, 15]. In these studies, the importance of inhomogeneous surface charge 
distributions has been recognized, and for a number of different charge-distribution models and 
parameters the ionic distribution functions as well as forces between charged surfaces have been 
calculated. In this paper we use a simple model for a charge-modulated surface which depends 
on a single geometric parameter and includes as limiting cases both smeared-out and delta- 
peaked charge distributions. We study the distribution of counterions without added salt at a 
single charged surface for all different values of the electrostatic coupling parameter (depending 
on temperature and counterion valence) and for different degrees of surface charge modulation, 
both using Monte-Carlo (MC) simulation methods and the recently introduced strong-coupling 
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(SC) theory [16, 17, 18]. Our study therefore encompasses different experimental situations such 
as mono/multivalent ions at surfaces with discrete charged chemical groups, as well as charged 
colloids at microscopically charge-modulated surfaces. 

As we will demonstrate, the assumption of smeared-out charges dramatically breaks down 
when the charge-modulation at the surface is pronounced, but also for moderate substrate- 
charge modulation when the electrostatic coupling is large (i.e. for highly charged surfaces, 
low temperatures, or multivalent counterions): in both cases the counterions become highly 
correlated with the surface charges and tend to form a two-dimensional, laterally ordered layer 
close to the surface. A convenient measure for the effects of substrate charge modulation is 
the counterion contact density, i.e., the laterally averaged counter-ion density at the substrate 
surface, since it is known exactly in the smeared-out case and can be easily determined from 
simulations [16]. We find that the contact density for charge- modulated substrates can be 
much larger than for the smeared-out case, which is in agreement with recent experimental 
measurements on highly charged surfactant monolayers [19]. Quite surprisingly, substrate 
charge modulation tends to have a more drastic effect on the counter-ion distribution than 
fluctuations and correlations, which have been the subject of numerous recent studies (see [16] 
and references therein). As our numerical results show, the SC theory describes the counterion 
distributions quantitatively in the SC limit and in particular in the limit of pronounced surface- 
charge modulation. Wc also consider a dielectric-constant jump at the substrate, as relevant 
for charged biological and colloidal surfaces[20, 21, 22, 23]. In contrast to the mean-field 
Poisson-Boltzmann (PB) approach, which severely fails even in the smeared-out case, our SC 
approach compares well with numerical results and helps to understand the intricate interplay 
of dielectric-jump and charge-modulation effects. 

The substrate charge distribution is modeled, for simplicity, by point charges distributed 
on a square lattice with lattice constant a. The point-like oppositely charged counterions are 
confined to the positive half-space with z > 0, the minimal distance between surface charges 
and counterions is given by D, as shown in Fig. la. The dielectric constant of the positive 
half-space e> is allowed to be different from e<. The system is globally neutral, and only 
the fixed surface charges and counterions are present, i.e., no salt is added. In principle, 
one could interpret D as the sum of the surface-ion and counterion hard-core radii. In our 
model, however, we use the ratio D/a more generally to control the degree of substrate-charge 
modulation, which allows to describe a variety of experimental situations and systems by 
a single parameter. The limit D/a — > oo is equivalent to a smeared-out surface-charge 
distribution, while D/a — > corresponds to a delta-peaked distribution. 

The typical height of the counterion layer in the smeared-out limit defines the Gouy- 
Chapman length fi = l/(2TrqlB&s), where a s = Q/a 2 is the number charge density at the 
wall, q and Q are the valences of counterions and surface ions, and £b = e 2 / ^Sye^ksT 
is the Bjerrum length (the distance at which two elementary charges interact with thermal 
energy, UbT). In the following, we will rescale all lengths by \i according to r = r/fi. For the 
simple double layer (smeared out charges at the surface and no dielectric jump) the coupling 
parameter S = q 2 (-B/ M = 2irq 3 £ B a s is the only parameter in the problem. In the limit S — > 
(weakly charged surfaces, low-valence counterions, or high temperatures) the PB theory is 
asymptotically exact, while in the opposite limit S — ► oo (highly charged surfaces, high- valence 
counterions, or low temperatures) the SC theory is exact [16]. The present model is in addition 
characterized by D/a, the dielectric-constant ratio A = (e> — £<)/(£> + e<) and the valence 
ratio q/Q between counterions and surface charges. 

The Coulomb interaction in the presence of dielectric discontinuities is given by the solution 
of the Poisson equation with the appropriate boundary conditions [24]. For a single dielectric 
jump located at z — as shown in Fig. la, the electrostatic energy reads 
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Fig. 1. - a) In our model a square lattice of fixed surface charges with lattice constant a is located in 
the xy-pl&ne at a distance D from the counterion containing half space, b) Simulation snapshots and 
normalized lateral density profiles p yz (x) (averaged in the j/z-plane) for a coupling constant S = 10, 
vanishing dielectric jump (A = 0) and valence symmetry between surface and counterions (Q = g) 
for D/a — 0.06 (left) and D/a — 0.24 (right), showing strong and weak lateral counterion ordering, 
respectively. In the snapshots, the surface ions are located beneath the lattice nodes. 
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The first sum corresponds to the interaction between pairs of counterions at positions 
(taking into account the so-called "image charges" through the term proportional to A), the 
second sum is the interaction between counterions and fixed surface charges (the sum over a 
corresponds to all surface- ion lattice vectors R Q ) and the third sum is the interaction between 
counterions and their "images." In our MC simulations we use Eq. (1) for a typical number 
of 100 counterions in conjunction with periodic boundary conditions (implemented by the 
Lekner-Sperb technique [2 5]) in order to minimize finite size effects[16]; the results reported 
here are always taken from runs with 10 6 Monte Carlo steps per particle. In Fig. lb we show 
the effects of varying surface-charge modulation: For a ratio D/a = 0.24 (to the right) the 
counterion snapshot shows a rather irregular configuration, with no or little visible correlation 
between surface charges (located beneath the nodes of the square lattice) and counter ions. 
For a ratio D/a = 0.06 (to the left), on the other hand, the counterions are strongly correlated 
with the square lattice of the surface ions (two counterions obviously have, driven by thermal 
fluctuations, escaped from their assigned lattice positions). This is reflected by the normalized 
lateral counterion distribution p yz (x), which for D/a — 0.24 oscillates weakly around its mean 
value of unity, while for D/a = 0.06 this distribution function is strongly peaked at the 
positions of the surface ions. 

Before we present the density profiles as a function of the distance from the substrate, we 
quickly recapitulate the SC theory. At leading order of an expansion in inverse powers of 
the coupling strength 5, which is equivalent to a virial expansion, the normalized counterion 
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Fig. 2. - Counterion density profile p xy (averaged in the xy-pla,ne) as a function of z = z/p for q = Q 
and without dielectric jump (A = 0). a) MC results (symbols) and SC predictions (full lines) for 
S = 10 and D/a = 0.06 (diamonds and inset), D/a — 0.12 (plus symbols), D/a = 0.24 (crosses). 
The circles, dashed and dotted lines denote, respectively, the MC data, PB and SC predictions for 
D/a — oo. b) MC results (symbols) and SC predictions (full lines) for D/a = 0.12 and H = 1 (stars), 
S = 10 (plus symbols) and E = 100 (squares and inset). Error bars are smaller than symbol size. 



density distribution is given by[18] 
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The factor A is the fugacity of counterions, determined by the normalization condition 
J dr p(r) / L 2 — 1, where L is the lateral system size, equivalent to the condition of global 
elcctroneutrality[18]. The function u, given by 
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is the electrostatic potential created at position r by the surface ions at lattice positions R Q 
(with Z a = —D) plus the counterion-image interaction. In previous studies we demonstrated 
that the SC approach becomes quantitatively accurate for coupling strengths larger than 
S w 100 [16]. The discreteness of surface charges tends to decouple different counterions from 
each other, as witnessed by the snapshots shown in Fig. lb. One would therefore expect 
corrections to the leading term of the systematic expansion in Eq.(2), which come from 
correlations between counterions, to be weakened and the SC approach to perform even better 
in the presence of modulated surface charges, as indeed borne out by our data. In the limit 
D/a — » oo, the smeared-out case, Eq. (3) reduces to u(r) = (1 + A)z + H A/(42). 

We first discuss our results in the absence of a dielectric jump (A = 0) and identical valences 
of surface and counterions (Q — q). In Fig. 2a we fix 3 = 10 while D/a is varied; in Fig. 2b 
we fix D/a — 0.12 while S is varied. Clearly, the counterion density profiles p xy (z) (which 
are averaged in the xy-plane) are very sensitive to both the coupling constant S and the 
ratio D/a. For 5 = 10 and D/a = 0.24 (crosses in Fig. 2a) the discretization has a small 
effect, as the data almost coincide with the MC results for the smeared-out case D/a — > oo 
(circles). With the rescaling of Fig. 2a, the difference between the smeared-out SC and PB 
profiles (dotted and broken lines, respectively) is in fact rather small compared to the effects 
of surface-charge modulation, and the smeared-out data (circles) are somewhat in between 
the SC and PB predictions, demonstrating that E = 10 is in the crossover regime between 
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Fig. 3. - Contour plots of constant contact counterion density p xy (z — 0) as a function of the coupling 
constant 3 and the ratio D/a for A = and Q = q. The lines follow from the SC theory via Eq. (2), 
while the open circles are obtained from MC data via extrapolation leading to p~xy(0) = 2.9, 5.9, 7.8, 
9.6 and 27.8. For D/a — * oo as well as B — » the smeared-out assumption for the surface charge 
distribution, characterized by pxy(0) = 1, is valid. 

strong and weak coupling[16]. For D/a — oo and A = the contact-value theorem predicts a 
rescaled contact density of unity, p xy (z = 0) = 1, as confirmed by simulations at various values 
of S[16]. However, as D/a becomes smaller there is a greater accumulation of counterions in 
the immediate vicinity of the charged surface, and the laterally averaged contact density can 
be several times larger than unity. This is reflected by a high lateral correlation between 
counterions and the surface charges, see Fig. lb. As D/a becomes smaller, the SC predictions 
as defined in Eqs. (2) and (3) (solid lines in Fig. 2a) show progressively better agreement with 
the MC results. Similarly, in Fig. 2b for fixed D/a — 0.12, the contact density rises while the 
agreement between MC data and SC predictions becomes quantitative as S is increased. 

These results are summarized in Fig. 3, where contour plots of constant contact density 
Pxy(z = 0) are shown as a function of S and D/a. The lines were obtained from the SC theory 
via Eq. (2) and agree well with the MC results (open circles), obtained from density profiles 
via extrapolation. In the limit D/a — > and finite coupling 3 the counterions collapse onto 
the surface charges and the density at contact diverges, while for fixed D/a and in the limit 
S — > oo the contact density saturates at a finite value. As D/a grows (or 2 becomes smaller), 
the smeared-out assumption becomes valid and p X y(0) approaches unity. 

We finally lift the restriction to symmetric surface/counterion valences (q = Q) and to 
vanishing dielectric jump (A = 0). In Fig. 4a we show density profiles for SQ 3 /g 3 = 10, 
D/a = 0.12, A = with varying valence ratios of q/Q = 1 (plus symbols, already featured in 
Fig. 2), q/Q = 2 (stars), q/Q — 4 (triangles) and q/Q = 16 (squares). The inhomogeneity of 
the surface charge distribution becomes more important as the counterion valence increases. 
The SC results (solid lines) capture this trend and quantitatively agree with simulation data. 
The contact density p xy (0) saturates at a constant value as q/Q — > oo (see inset). For the 
case where a dielectric jump is present, we choose £> = 80 and e< = 2, leading to A = 0.95 
and corresponding to a charged hydrocarbon substrate in contact with water. Since A is 
positive, the interaction between the counterions and their images is repulsive and competes 
with the attractive interaction due to the surface charges, see Eq. (1). Consequently, the 
single-ion potential u(r), Eq.(3), exhibits a minimum at a finite distance from the wall, and 
one expects the maximum in the density profile to be displaced from the wall, in agreement with 
previous results[20, 21] (however, one should note that the resulting density profile depends 
strongly on details of the geometry and position of the dielectric boundary). Fig. 4b shows 
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Fig. 4. - a) Counterion density profiles p xy (z) for HQ 3 /q s — 10, D/a = 0.12, A = and different 
values of the valence ratio q/Q. Full lines correspond to the SC predictions and symbols to MC results 
for q/Q = 1 (plus symbols), q/Q — 2 (stars), q/Q = 4 (triangles) and q/Q — 16 (squares), b) MC 
results (symbols) and SC predictions (full lines) for the counterion density profile for S = 10, q — Q 
in the presence of a dielectric jump, A = 0.95, for D/a — 0.06 (diamonds and inset), D/a = 0.12 
(plus symbols) and D/a = 0.24 (crosses). The circles, dashed and dotted lines denote, respectively, 
MC data, PB and SC predictions in the smeared-out limit, D/a — » oo. 



the counterion density profile for S = 10, q = Q, A = 0.95 and the same values of D/a as 
in Fig. 2a, namely D/a = 0.06, 0.12 and 0.24. The corresponding SC predictions (shown as 
solid lines) locate correctly the position of the maximum in the density, though the agreement 
is not as good as without the dielectric jump for the larger values of D/a. The data for the 
smeared-out case, D/a = oo (open circles) exhibit a maximum close to z = v/HA/(4(l + A)), 
the SC prediction, although the SC profile for this case (dotted line) overestimates the amount 
of counterions close to the wall. Obviously, the SC theory does not take into account the 
repulsion felt by a counterion due to the presence of the images of other counterions, which 
explains why SC theory performs at the same value of 2 better in the absence of a dielectric 
jump (compare Fig. 2). Higher order corrections to the asymptotic SC theory would therefore 
be needed to quantitatively match the data. As is well-known, PB theory (shown as a broken 
line) is insensitive to the presence of a dielectric jump in the smeared-out limit D/a — > oo and 
completely misses the shape of the density profile [20, 21]. 

To conclude, we have studied both analytically and numerically the distribution of coun- 
terions close to an inhomogeneously charged surface. The charge inhomogeneity affects the 
counterion distribution strongly for pronounced charge modulation as well as when the coupling 
constant 3 is large. As summarized in Fig. 3 in the absence of a dielectric discontinuity (A = 0), 
the laterally averaged counterion contact density is much larger than the smeared-out value 
at decreasing D/a (equivalent to a strongly charge-modulated substrate) and at increasing 
coupling constant S, and counterions become strongly correlated with the surface charges. 
The same effect is reflected in the counterion density profiles, which under such conditions 
decay faster to zero as one moves away from the substrate, see Fig. 2. 

This has direct experimental consequences: Very recently the counterion density profile 
(measured using ellipsometry) at a surfactant monolayer with a high surface charge density 
(determined via second-harmonic generation) could be described by PB profiles only by as- 
suming a fraction of counterions to be bound to the monolayer [19], similar to what we find. 
Our results suggest to reconsider the traditional picture of a Stern layer, where the counterion 
surface concentration is subject to electrostatic and specific surface attraction and packing con- 
straints, but where the influence of lateral surface charge modulation (which is always present 
experimentally and leads to sizable effects as shown here) is typically neglected. One should 
emphasize the experimental relevance of our parameters: for a system with a s ~ 1/257 A -2 
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and divalent counterions in water (which at room temperature leads to S = 10), the surface 
ions are on average at a distance a ~ 16 A from each other. For a minimal distance between 
surface ions and counterions of D ~ 2 A (noting that at close contact the hydrated water shell is 
stripped off) the ratio D/a = 0.12 is obtained; as our results demonstrate, the inhomogeneous 
character of the surface charges leads to pronounced deviations from the smeared-out case. 

Other relevant factors, which we will consider in future work, include the possibility of the 
counterions to actually penetrate between the surface ions, different surface charge patterns 
(annealed and various quenched geometries), and additional interactions between surface 
groups and counterions (like van der Waals or solvation forces) and between counterions 
themselves (most importantly hard core interactions, which to leading order within the SC ap- 
proach simply shift the density profile away from the surface). The advantage of the currently 
employed model, though being quite crude and neglecting a number of important effects, is 
that its simplicity allows for a global analysis, encompassing different and limiting values of the 
coupling strength and the degree of substrate charge modulation. In a complementary study, 
similar numerical calculations have been recently analyzed within a modified PB approach [26]. 

We thank H. Motschmann and R. Tcppncr for useful discussions. AGM acknowledges 
financial support from the DFG Schwerpunkt Polyelektrolyte. 
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